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Abstract. We introduce the most general version of Dubrovin-type equations for divisors on a 
hyperelliptic curve fC B of arbitrary genus g 6 N, and provide a new argument for linearizing the 
corresponding completely integrable flows. Detailed applications to completely integrable sys- 
tems, including the KdV, AKNS, Toda, and the combined sine-Gordon and mKdV hierarchies, 
are made. These investigations uncover a new principle for 1 + 1-dimensional integrable soliton 
equations in the sense that the Dubrovin equations, combined with appropriate trace formulas, 
encode all hierarchies of soliton equations associated with hyperelliptic curves. In other words, 
completely integable hierarchies of soliton equations determine Dubrovin equations and asso- 
ciated trace formulas and, vice versa, Dubrovin-type equations combined with trace formulas 
permit the construction of hierarchies of soliton equations. 



1. Introduction 

The purpose of this paper is to re-examine Dubrovin equations for divisors on hyperelliptic 
Riemann surfaces and to underscore their exceptional role in connection with completely integrable 
hierarchies of soliton equations. 

Starting from four representative hierarchies, the Korteweg-de Vries (KdV), Ablowitz-Kaup- 
Newell-Segur (AKNS), Toda lattice (Tl), and the combined sine-Gordon and mKdV (sGmKdV) 
hierarchy we derive a new argument for linearizing the corresponding completely integrable flows. 
As a result of these investigations we show that a proper combination of Dubrovin equations and 
trace formulas involving auxiliary divisors on hyperelliptic curves encodes all information on the 
underlying completely integrable hierarchy of soliton equations. 

In Section || we briefly review some basic facts on hyperelliptic curves and establish the no- 
tation used throughout this paper. Section |^ provides a "crash course" into the four different 
hierarchies closely following the detailed accounts in (ll), Q, f35j , and In particu- 

lar, we outline an elementary polynomial, recursive approach to these hierarchies, as originally 
developed by S. I. Al'ber §, @ in the KdV context, and introduce the corresponding divisors 
on /C g associated with their algebro-geometric solutions. Section [| relates the polynomial recur- 
sion relation approach with elementary symmetric functions (and functions derived from them) 
of (A\(x, t r ), . . . ,/j, g (x,t r ), where jij (x, t r ) are certain analogs of Dirichlet-type eigenvalues of the 
corresponding Lax operator. In Section we study Dubrovin equations and, based on the results 
of Section [|, provide a new proof of the fundamental fact that a change of coordinates effected by 
the Abel map straightens out the flows on the Jacobi variety J{TC S ) of K s . In Section Jq we briefly 
recall the theta function representations of algebro-geometric solutions and Section u\ illustrates 
our results by deriving interesting connections between the KdV fl and sG equations, and AKNS 
and Toda hierarchies, respectively. These connections establish the fundamental role played by 
the Dubrovin equations as the common underlying principle for hierarchies of soliton equations. 
In particular, our formalism establishes an isomorphism between the class of algebro-geometric so- 
lutions of these integrable systems. Finally, Appendix |A| collects some useful results in connection 
with Lagrange interpolation formulas. 
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We emphasize that our results are not necessarily restricted to hyperelliptic curves. In partic- 
ular, the approach of this paper applies to Boussinesq-type curves using the polynomial recursion 
formalism for the Boussinesq hierarchy developed in [ fl9| , pp| . 

Depending perhaps a bit on one's taste, the results of this paper may at first sight appear 
somewhat discouraging as they clearly shift the emphasis from individual hierarchies of soliton 
equations toward Dubrovin-type equations. On the positive side, however, they establish the 
Dubrovin equations as a universal object underlying all hierarchies. 

We note that our approach to completely integrable soliton equations is close in spirit to that 
developed by M. S. Al'ber and S. I. Al'ber in a series of papers (see, e.g., and the references 

therein) . While their approach focuses on algebraically integrable systems and hence on a Hamil- 
tonian formalism with associated action and angle variables, our approach concentrates on how 
a combination of elementary symmetric functions of Hi(x, t\), . . . , fj, g (x, t r ) and certain trace for- 
mulas generate completely integrable hierarchies of soliton equations and their algebro-geometric 
solutions. 

A different series of papers closely related to our investigations and focusing on characterizing 
real-valued solutions of various soliton equations was published by Chen, Chin, Lee, Neil, Ting, 
and Tracy |5^|-[^l), |63|-[[35). These authors, however, appear to be unaware of the prior work of 
M. S. Al'ber and S. J. Al'ber in this field. 

Finally, we stress that the use of elementary symmetric functions and hence of trace formulas 
in terms of Dirichlct eigenvalues has a long history in the context of integrable equations. In fact, 
as early as 1975, Flaschka |30) characterized the real-valued periodic potentials q with finitely 
many stability intervals of the associated Schrodinger operator —d 2 /dx 2 + q in L 2 (R;dx) as sta- 
tionary solutions of the KdV hierarchy using (regularized) trace relations for Dirichlet eigenvalues 
associated with q and the underlying periodicity interval. 



2. Hyperelliptic curves 
Fix N £ No- We briefly review hyperelliptic Riemann surfaces of the type, 

N 

T N (z,y) = y 2 - R N+ i(z) = 0, R N+1 (z) = _Q (z - E m ), 

m=0 

{-E m } m=0 ,...,JV C C, E. m ^ E m , for m ^ m'. (2.1) 



The material of this section is standard and can be found, for instance, in |29]|. The curve (2.1) is 
compactified by adding one point at infinity if N is even, and two points P (x>+ and Poo_ if N 
is odd. 

One introduces an appropriate set ofQ [N + lj /2 nonintersecting cuts Cj joining E m ^ and 
E m i(j) and joining En and oo if N is even. Denote 

C= |J C h CjnC k = ®, j^k, (2.2) 
je./u{oo} 

where J C {1, . . . , [N + 1J/2}. Define the cut plane 

n = C\C, (2.3) 

and introduce the holomorphic function 

/ JV \ V2 

i? w+1 (-) 1/2 : n^C, z^[Y[{z-E m )\ (2.4) 



on II with an appropriate choice of the square root branch in (2.4). Define 

'{Poo} i 

{Poo+,Poo_} for N odd 



M s = {{z,aR N+1 {zf' 2 )\z^C, ae{±l}}u|^°° } „ 1 ^ ™' (2.5) 



[x\ = sup{j/ e Z | y < x}. 
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by extending _Rjv + i(-) 1 / 2 to C. The hyperelliptic curve K. g is then the set A4 g with its natural 
complex structure obtained upon gluing the two sheets of M g crosswise along the cuts. Finite 
points P on JC g are denoted by P = (z, y), where y(P) denotes the meromorphic function on K, g 
satisfying Tn{z, y) = y 2 — Rn+i{z) = 0; K, s has genus g — [N + lJ/2. 

A basis of g linearly independent holomorphic differentials on lC g is given by z l ~ x dz/y{P) for 
£ = 1, . . . , 0, and we introduce 

duj( p ) = *51 c i /, l( p\ ' 3 = l,---,fl, (2-6) 



: //(/') 



with normalization, 



/ duj=8 jlk , (2.7) 

Ja k 



where {dj, bj}j =1 is a homology basis for JC g . 
Define the matrix r = (jj.t) by 



Tj,e = [ duj e . (2.8) 



Then Im(r) > and t^j — nj. 

The Riemann theta function associated with /C and the given homology basis {fflj, bj}jLi> by 
definition, reads 

0{z) = ^2 cxp (27ri(n,z) +m(n,Tn)), zG<C fl . (2.9) 

nGZ8 

We fix a base point Po on /C fl and define the Abel map A Pq by 

A Pa (P) = ( [ du) U ...,( dcu g ) (modL fl ), P G JC g , (2.10) 

J Pa JPa 

with period lattice 

L g = {n + rm | n,m G Z 8 }. (2.11) 

Similarly, we introduce 

a Po : Div(/C B ) - J(/C fl ), D ^ 2Po (P) = ^ P(F)4(P), (2.12) 

where Div(/C ) and J(fC g ) = C B /L g denote the set of divisors on /C and the Jacobi variety of /C fl , 
respectively. 

In connection with divisors on )C g we shall employ the following (additive) notation, 
Vq q = V Qo +Vq, Vq = V Qi +--- + V Qn , etc., 

for Q = (Qi,...,Q„) G <J n JC g , (2.13) 



where for any Q G /C fl , 



p »w> = J £^£\ W}l (2 ' 14 > 



and <J n K, g denotes the nth symmetric product of 1C g . 
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3. The hierarchies 

We give a brief presentation of the KdV, AKNS, sGmKdV, and Toda hierarchies based on a 
polynomial, recursive approach. The material of this section originated with work of S. I. Al'ber 
@, U (see also §-§). It has been further developed in glj , || , [|| , |E) , and [§§ and we closely 
follow the latter sources. Common to all these hierarchies is that one can naturally associate with 
each one of them a hyperelliptic curve as described in the previous section. 

The KdV hierarchy. The Lax pair consists of a second-order linear differential expression L 
of Schrodinger-type, 

d 2 

L (h) = —^ + V{x,t s ), (3-1) 

and a differential expression P2 B +i(t g ) of order 2g+ 1 defined recursively as follows. Let {/j}jeNo 
be given by 

/o = l, /,, \.f, • \ f, :., I \v x fj I- jGN. (3.2) 

Explicitly (fj = fj(x,t B ), (x,t B ) 6 R 2 ), 

/i = |V + ci, h = -\V xx + lV 2 + Cl \V + c 2l etc., (3.3) 
where c, € C are integration constants. Then one defines 



5 d 1 



_\a-j 

' dx 



G No, (3.4) 



and using the definition of fj in (3.2), one finds that the commutator of L(t g ) and P2 g +i(t g ) is in 
fact a multiplication operator. Indeed, the Lax commutator representation reads 

M*fl) ~ [ p 2 fl+ i(i ),i(« B )] = ^ ~ 2f g+hx (x,t g ) = KdV g (V) = 0. (3.5) 

Explicitly, one obtains for the first few KdV equations, 

KdV o (V) = V to -V x = 0, 

KdVi(K) = K a + \V XXX - \VV X - c x V x = 0, (3.6) 

Kdv 2 (y) = v t2 - ^v xxxxx + \vv xxx + \v x v xx 

- fV 2 V x - c 2 V x + cx{\V xxx - \W X ) = 0, etc., 

where, of course, KdVi(V) is the KdV equation. 
Next define the polynomial^] of degree g in z 

B 3 

F B (z,x,t g ) = ^/ fl _j(a;,i fl )z J = - Hj(x, t )) (3.7) 

3=0 3=1 

implying 

-2V ta = F g . xxx - 4(1/ - z)F B . x - 2V x F g . (3.8) 
In the special stationary case, defined by Vt — 0, this integrates to 

\f b , xx F b - - A Fl x -(V- z)F 2 = R 2B+ i(z). (3.9) 
Here i?2 e +i is a monic polynomial of degree 2g + 1 with zeros {Eq, . . . , E2 B }. Hence, 

R2 5 +i{z) = [I ( z {^}™=o,..,ai c C - ( 3 - 10 ) 



2 The zeros fj,j(x,t s ) of -F B turn out to be eigenvalues associated with L(t B ) and a Dirichlet boundary condition 
at the point 
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The hyperelliptic curve IC g is defined in terms of the stationary KdV hierarchy obtained by 
considering a ^-independent function V = V{x), resulting in 

[P 2fl+1 ,L] = 2/ B+M = 0. (3.11) 



The classical Burchnall-Chaundy theorem 1 12 , |13| , 1 14 (see also [|56j , |66| ) states that commuting 
differential operators are algebraically related. In the present context one finds 

Pi B+ i=R2 e+ i(L) (3.12) 

and thus the hyperelliptic curve K, g of genus g is of the type y 2 = i?2 B +i(z), with N = 2g even 
when compared to Section ||. 

For later purpose we quote the following asymptotic high-energy expansion^ (see ||l]], |3(|, (57j) 



^J3 i Z > X > ) 

2R 2g +i(z) 1 / 2 z^oo 2zV 



1/2 fj( X ' t$) 2 



(3.13) 



where fj denotes the homogeneous coefficients fj with all integration constants equal to zero, 
ct = 0, 1 > 1, that is, 



/ = /0 = l, /j=/i| C f=0, £=l,...,j, JGN. 

We also introduce the following fundamental meromorphic function on K. g , 
iy{P) + \F^x{ z ^ x ^s) -Hg+i{z,x,t g ) 



<f>(P,x,t B ) 



F B (z,x,t B ) iy(P) - ±F BjX (z,x,t B )' ^ x,t ^ E 



(3.14) 



(3.15) 



(the second equality in ( 3.15| ) serving as a definition of the polynomial H g+ i of degree g + 1 with 
respect to z) and the time-dependent Baker-Akhiezer function ^p{P, x, xq, t g , io,g) on JC g \ {Poo}, 

ip(P,x,x 0> t B ,t 0<B ) 



exp 



ds (F B (z,x ,s)4>(P,xo,s) - lF BtX (z,x ,s)) + / dx' 4>(P,x' ,t B ) 

(x,x ,t B ,t 0tB ) e M 4 . 



The divisor (<p(P,x,t B )) of 4>(P,x,t B ) is given by 

(<f)(P,X,t B )) ='Dp (x,t si )v{x,t s ) - r Dp oa i i {x,t s )i 

where 



(3.16) 
(3.17) 



Ai(M B ) = (fi j (x,t B ),-^F BtX (n j (x,t B ),x,t B )), j = l,...,g (3.18) 
denote the Dirichlet divisors! and 



^(MjO^OMb) = {vo{x,t B ), ■ ■ ■ ,0 s (x,t B )) E a e+ /C B , 

h(x,t B ) = {y t [x,t 9 ), |P , x (^(a;,t B ),a;,i fl )), £ = 0, ...,g (3.19) 
abbreviate the Neumann divisors^] derived from the zeros of -ff B+ i(z, x,t g ), 



H g+1 (z,x,t g ) = Y[(z - v t {x,t B )). 



(3.20) 



The importance of <p in connection with divisors on hyperelliptic curves was recognized by Jacobi 
[ 41 1 and applied to the KdV case by Mumford Q, Sect. IIIa.1 and McKean [|9| (see also pi) , 
|56|). 



3 For an appropriate choice of the sign of the square roots in (3.13 ), the left-hand side of (3.1c ) equals the Green's 
function on the diagonal (i.e., for x = x 1 ) of L(t s ). 

4 The zeros Vf(x, t g ) of -ffg+i turn out to be eigenvalues associated with L(t g ) and a Neumann boundary condition 
at the point 
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The AKNS hierarchy. The Lax pair consists of a Dirac-type matrix-valued differential 
expression 

M(t - )=i C(4) ~ q -t } )' t>€R - <3 ' 21) 

and a matrix-valued differential operator Q B +i(t g ) of order g + 1. To define Q B +i(i B ) we proceed 
as follows. Define {f e (x,t 3 )} eeNo , {g e (x,t 3 )} eeNo , and {h e (x,t 3 )} eeNo recursively by ((x,t 3 ) G M 2 ), 

fo = -iq, 9o = l, h = ip, (3.22) 

i i 

fe+i = -=jt,x - iqgt+i, ge+i,x = pft + qhi, h e+1 = --h ttX + ipgt+i, £eN . 



^w-e( 5riV /rfriW^ ^no, /- 1 =fc_ 1 =o, (3.23) 



The 2x2 matrix Q fl +i(t B ) is then defined by 

8+1 

-K-eih) 9a+i-t(t g ) y 
and one verifies that the commutator of Q s +i(t g ) and M(t g ) becomes 

Consequently, the Lax commutator representation for the AKNS hierarchy reads 

- », +1 ( t .),M«,)] . - . o, 

g G N . (3.25) 

The first few equations equal, 

AKNS (p,g) = f^- fc + Cl( - 2 t ) )=0, 

V ?t D -fe + ci(2i 9 ) y 

/ Pti + Ip** - *P 2 <? + c i (~Px) + c 2(-2ip)\ 
AKNSi(p,g)= 2 =0, etc. (3.26) 

\ Qh ~ \q_xx + ipq 2 + ci (-£fe) + c 2 (2iq) J 

Next, define polynomials F B , Gj+i, and H s with respect to z G C, 

s 9 
F 8 (z,a;,t B ) = ^/ fl _*(x,t B )^ = t„) JJ(^ - /^(a;, *„)), 

s+i 

G Q+ i{z,x,t g ) = ^2g g+1 - e (x,t s )z e , (3.27) 

8 8 

if B (z, x,i 8 ) = ^2h g ^ e (x,t g )z e = ip(x,t g ) Y[(z - fj(x,t B )). 
e=o 3=1 

In the special stationary case, where pt e = <Zt e = 0, we infer from the recursion ( 3.22| ) that 

(G 2 3+1 -F s H g ) x =0 (3.28) 

and hence 

G B + 1 - F s H a = R2 9+ 2(z), (3.29) 
where R 2b +2 is a monic polynomial of degree 2g + 2 with zeros {Eq, . . . , i?2 B +i}. Thus, 

2fl+l 

R2g+2(z) = ]J ( Z ~ {^m}m=0,...,2fl+l C C. (3.30) 

7TI — 
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The stationary case determines the hyperelliptic curve JC g of genus g of the type y 2 = R2 B +i(z), 
with N = 2g + 1 when compared to Section |^. If p = p(a^) and g = g(a;) are stationary solutions 
of the AKNS equation, 

[Q B+ i,M] = 0, that is, f g+1 = h B+1 = 0, (3.31) 
Burchnall-Chaundy's theorem implies that 

Q 2 B+1 + R 2S+2 (M) = 0. (3.32) 
By studying the Green's matrix of M one finds the following asymptotic high-energy expansion 

n, 

f'^f Wa = -f>(M fl K fe , (3.33) 

i?2 B+2 (z) 1/2 2 ^ 

H g (z, X, t g ) 1 ^ t / . \ —k to OA\ 

p = -Z^ h k{x,t B )z (3.34 

for an appropriate determination of the square roots in (3.33) and (3.34). Here fj and hj denote 
the homogeneous quantities with vanishing integration constants q, £ > 1, that is, 

fo = h = -iq, fj = /iU=o, <=1 AjeN, (3.35) 

h = h = ip, hj = hj\ ct =o, £ = l,...,j, j €N. (3.36) 

We also record the functions 

, (p , n +gs+l^l^j) = ~ h b( z ^ x ^b) 

n ' ' bJ M b ) t,(P)-G B+ i(z,Mg)' 

P=(z,y)e/C B (3.37) 

and the Baker- Akhiezer vector, 



( ih(P,x,x ,t B ,to, g ) 

\i) 2 (P,x,XQ,t B ,t , B ) 

ipi(P,x,x 0l t g ,t ^ g ) = cxp I / da;' (-iz + q(x' ,t B )<f>(P,x' ,t B )) 

\J Xq 



x 

+i / ds(F B (z,x ,s)(f>(P, x ,s) - G B+ i(z,x ,s)) J , 

^(-P^^o^g^o.s) = <t>{P,x,t g )ipi(P,x,x ,t g ,t , B ), (3.39) 

The divisor of <fr(P, x,t g ) is given by 

{(j){P,X,t B )) = T> Poo+ p (Xits) -T>p oa _p,( x ,t s ), ( 3 - 40 ) 

where 

AOMa) = ■ ■ -jMsOMb)) G cr fl /C B , 

p,j(x,t g ) = (fj, j (x,t B ),G g+1 (n j (x,t B ),x,t B )), j = l,...,Q (3.41) 

and 

^(^ifl) = ("j^i^i-Gj+i^j^ij),!,^)), j = l, ...,fl. (3.42) 

The Toda hierarchy. Let (S ± f)(n) — f (n) = f(n ± 1), n 6 Z denote the shift operation 

on the lattice Z. The equations in the Toda hierarchy are continuous in time and discrete in the 
space variable. The Lax pair consists of the second-order difference operator 

L(t g ) = a(t g )S+ +a-{t g )S- -b(t s ), t g <ER (3.43) 
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and a difference operator P 2g+2 (i g ) of order 2g + 1 

P2 B+2 (t B ) = -L(t 5 )° +1 + J2(9j(t B ) + 2a{t 5 )f 3 (t,)S + )L{t 5 y-i + f B+ i(t 3 ), 

j=o 

g G N . (3.44) 

Here {fj(n, t B )}j £ n ancl {9j( n ^ B )}jefh satisfy the recursion relations, 
fa = 1j 9o = -ci, 

2/i+i + <?j + 3/ + 26/, =0, je No, (3.45) 
fli+i - 57+1 + 2 (« 2 /+ - (a-) 2 f-) + b(g 3 - gj) = 0, je N . 
The Lax commutator representation of the Toda hierarchy then reads 

it, - [P2 B+ 2(t s ),L(t 3 )} = T1 B (a, 6)rS+ - T1 B (a, 6) 2 + Tl g (a~, b~)iS~ = 0, 

g e No, (3.46) 

where 

Tl g (a, 6)x = a to + a(. 9 + + , 9g + /+ +1 + f g+1 + 26+/+) = 0, 
Tl (a, 6) 2 = 6 to + 2(6(. 9g + / g+1 ) + a 2 /+ - (a-) 2 f B + b 2 f g ) = 0. 
This is equivalent to 

Tl B (o,6) = (Tl fl (a,6) 1 ,Tl B (a,6) 2 ) t = 0, g e N . (3.48) 
The first few equations of the Toda hierarchy equal, 

/a t2 -a((a+) 2 -(a-f + (b+f-b 2 )\ ( -a(b-b+) 
U ' ' \b t2 - 2a 2 (b+ + 6) + 2(a-) 2 (6 + b~) J 1 ^-2 ((a") 2 - a 2 ) 
etc. 

Next, define 

8 

F s {z,n,t g ) =^r f z j f - j (n,t ) = Y[(z - fJ,j(n,t B )), (3.50) 

3=0 3=1 
3 

G B+1 (z,n,t g ) = -z 3+1 + '£z 3 g a - j (n,t g ) + f g+1 (n,t g ). (3.51) 

3=0 



(3.47) 



In the special stationary case, defined by a tfl = b tg = 0, the recursion formulas (3.45) then imply 

G B+1 - 4a 2 F B F+ = G B+1 - 4( a -) 2 F g F 3 = i? 2fl+2 (z), (3.52) 

where i? 2g + 2 (z) is a lattice constant. By inspection, i? 2g+2 (z) is a polynomial in z of degree 2g + 2 
with zeros {Eq, . . . , Z? 2g+ i}, that is, 

20 + 1 

R2 B +2(z) = Y\_ ( z ^m), {Sm}m=0,...,2 8 +1 C C. (3.53) 
m— 

Consider now the stationary hierarchy where a = a(n) and b = 6(n) satisfy [P20+25 = 0, or 

Ka, + /,+■) + » 2 /, + - (OV, - + = o. 

Burchnall-Chaundy's theorem then states that 



% +3 = i? 2s+2 (L). (3.55) 
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Hence the hyperelliptic curve /C B of genus g reads y 2 = B,2 3 +2(z), and thus N = 2g + 1 is odd in 
the terminology of Section 0. 

Studying the diagonal Green's function of L(t g ) yields the high-energy expansion JTlj ] 



F {z,n,t B ) 



(3.56) 



with an appropriate choice of the radical in ( |3.56 ). Here fj and similarly gj, denote the homoge- 
neous coefficients fj and gj with vanishing integration constants ct = 0, I > 1, that is, 



/o — 1, fj — fj\c e =Q, 



■J, 3 e N, 



5o = 0, gj = g 3 \ Ce =o, £ = l,...,j + 1, j €N. 
Furthermore, define the meromorphic function 4>{P, n, t B ) on /C B by 



(f>(P,n,t B ) = 



-G s+1 (z,n,t g ) +y(P) _ -2a(n,t g )F B (z,n + l,t g ) 



(3.57) 
(3.58) 



P=(z,y) (3.59) 



2a(n,t 3 )F g (z,n,t g ) G g+ i(z,n,t B ) + y(P) 

using relation (|||). With the help of cj)(P,n,t g ) we define another meromorphic function 
tp(P, n, no, t g , to, g) on IC g , the Baker-Akhiezer function, by 



ip(P,n,n ,t g) t 0tg ) 



exp 



ds(2a(n , s)F fl (2, n , s)4>(P, n ,s) + G g+ i(z, n , s)) J x 

llm^no <^ P ' TO ' *b ) for n > n + 1, 
x < 1 for n = no, (3.60) 

.n^W,"*,*,)- 1 forn<n -l. 



The divisor (4>(P,n,t g )) of (j>(P,n,t g ) is given by 

(4>(P,n,t B )) = % + p(„+i, ff ) ~ 2? Poo_A(n,tg)> 

where 

£(m,* B ) = (Ai(m,t g ), . . . ,£ 8 (m,i B )) £ cr fl /C B , m e Z, 
fi 3 (n,t B ) = {fJ,j(n,t B ), -G B+1 (nj(n,t B ),n,t B )), j = !,...,$, 
fij{n + l,t B ) = (fij(n + l,t B ),G e+ i(pj(n+ l,t B ),n,t g )), j = 1, . 



(3.61) 

(3.62a) 
(3.62b) 
(3.62c) 



The sGmKdV hierarchy. The combined sine-Gordon and mKdV hierarchy is defined in 
terms of a zero curvature formalism as follows. Introduce the 2x2 matrices 



and 



J7 (z, a;, t B ) = — i 



. (\u x {x,t B ) 1 

■^u x {x, t B ) 



2 "<a;v~! "O; 
z 



(x,t g ) 6 



(a;,i fl )eIR 2 , GN O . 



H B (z,x,t g ) G B -i(z, x,t B ) y 
Then the zero curvature relation reads 

^-^ + [[7,^1=0, geNo, 

resulting in the equations 

u x t g (x,t g ) = -2iG B -i iX (x,t B ) - 2(H g (x,t g ) - F g (x,t g )), 
F B , x {x,t g ) = -iu x (x,t B )F g (x,t g ) - 2izG g - 1 (x,t B ), 
H BtX (x,t g ) = iu x (x,t g )H B (x,t B ) + 2izG g - 1 (x 1 t g ). 



(3.63) 

(3.64) 

(3.65) 

(3.66a) 
(3.66b) 
(3.66c) 
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Making the following polynomial ansatz 

S 8 

F B (z,x,t g ) = ^2f g _ j (x,t g )z J = Y[{z - Hj{x,t a )), (3.67a) 

3=0 3=1 
S B 

H s (z,x,t B ) =^2h g - j (x,t s )z :l = Y[(z - v>j(x,t B )), (3.67b) 

3=0 3=1 

0-1 

G-i{z,x,t g ) = 0, G s -!(z,x,t g ) = ^2g g ^ 1 ^ :j (x,t B )z 1 , (3.67c) 

3=0 

one concludes^ (see Q for a detailed discussion) (fj = fj(x,t g ), etc.) 

f g = ae- m , h g = (5e lu , a, (3 G C, g G N . (3.68a) 

/o = 1, ./'....- |/j • "'■// i . .• • "••.,/, :• j = 1, ■ ■ ■ ,0, G N, (3.68b) 

h = l, h j>x = ~hj-x, xxx + w-hj-i, x + ^ui- iX hj-i, j = 1, . . . ,g, g G N, (3.68c) 



where 



and 



«>± = ± Ziuxx), (3.69) 

3-1=0, gj = -(fj, x + iu x fj) = -(-h jtX +iu x hj), j = 0, . . . ,g - 1, g G N. (3.70) 



Explicitly, 

.go = -jWx, 3i = ti"^ + \uxxx - %u x , etc. (3-71) 

We also list a few coefficients in the homogeneous case where all integration constants eg,, I > 1 
vanish, 

/o=/o = l, 3o = 3o = -\u x , h =h = l, g £ N, 

h = fj\ Ce =o> 9j = 9j\ Ce =o> h i = ft i| C| =o' ( 3J2 ) 

^ = 1, ■ • ■ ,3, 3 = 1) • • • ,0 - 1, > 2, 
f g = f B = ae-™, h B = h B =0e iu , g G N . (3.73) 

Remark 3.1. The recursion for the sGmKdV hierarchy is anomalous compared to the other hierar- 
chies studied in this paper. For the KdV, AKNS as well as the Tl hierarchies the functions fj (and 
gj and hj where applicable) are defined by the same recursion formula for all j G N irrespective 
of the given genus g. However, for the sG hierarchy f g and h g are always given by (3.68a) for 
G No- This raises a compatibility problem in the recursion formalism. A proof of the solvability 
of the recursion can be found in 1 33 , Appendix C. 

The gth sGmKdV equation is then defined by 

B GmKdV 9 (tt(a;, t g )) = u x t B {x, t g ) + 2ig g ^ x (x, t g ) + 2( / 3e m ^ t «> - ae~ iu M) = 0, 

0GN o . (3.74) 

Explicitly, the first few equations read 

sGmKdVo(u) = u xto + 2{f3e m - ae~ lu ) = 0, 

sGmKdVi(w) = u xtl - iu xx + 2((3e tu - ae" 1 ") = 0, (3.75) 
sGmKdV 2 (w) = u x t 2 + i(u 3 x + 2u xxx ) x - ciiu xx + 2(f3e m - ae~ lu ) = 0, etc. 



J Observe that the recursions (3.68b) and (3.68c) arc iden tical with the KdV recursion replacing V by w±. Hence 
explicit formulas for the first few can be read off from ( h.g| ) upon changing V into w±. 
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Observe that a = = i/A and g = yields the well-known sine-Gordon equation in light- 
cone coordinates. Appropriate choices of a and (3 in (3.74) include the sinh-Gordon hierarchy, 
the corresponding elliptic equations, the Liouville model, as well as the modified KdV hierarchy 
(taking a = f3 = 0). 

In the stationary case, where u x t a — 0, we find 

■^(zG g -i(z,x) 2 + F g (z,x)H g (z,x) ) =0 (3.76) 



dx 

and hence 

zG^z^) 2 + F s (z,x)H g (z,x) = P 2g {z), (3.77) 



where P 2g {z) is x- independent. It is more convenient to define R2 S +i(z) = zP2 g {z) so that (3.77) 
becomes 

z 2 G B _i(z, x) 2 + zF g {z,x)H g (z,x) = R 2g+ i(z), (3.78) 
where R 2g +\ is a monic polynomial in z of degree 2g + 1 of the form 

fl2 fl +i(z) = l[(z-E m ), E = 0, E 1 ,...,E 2g eC. (3.79) 

m— 

This polynomial defines the hyperelliptic curve ICn of genus g by the relation y 2 — R2 g +\(z) = 
and hence N — 2g in the terminology of Section g. K, g is compactified by adding a point . 

Remark 3.2. Observe that the sGmKdV-curve is a special case of the KdV curve with the addi- 
tional constraint Eq = 0. 

Remark 3.3. In the stationary case the choice of a and (3 is constrained by the relation 

2fl 

a/3 = Y[E :j , (3.80) 

i=i 

as can be seen by inserting z = in ( |3.78 ), using ( 3.68a ) and ( |3.7S| ). 

We now return to the time-dependent formalism. Let (j){P,x,t g ) be the meromorphic function 
on /C B given by 

,/p , \ y{ p ) - zG a -i.{z,x,t a ) _ zH g {z,x,t g ) 

(x,t B ) eM 2 , P= (z, y ) e£ fl \ {*»«,}. 

Hence the divisor (cj)(P,x,t g )) of (f>(P,x,t g ) reads 

(<p(P,x,t g )) =^Q 0M (x,t a ) -^F^Cx.ts)' ( 3 - 82 ) 

with 

Kx,t g ) = (/ii(x,i B ), . . . ,/} B (a;,ig)) £ cr 8 /C B , (3.83a) 
fi 3 (x,t r ) = {nj(x,t r ),-^{x,t r )G g -i{tJ,j{x,t r ),x,tr)) £ K g , j = l,...,g, 
£(^)*s) = (AtOMs),---,^^,^)) Gcr fl /C B , (3.83b) 
i>j{x,t r ) = (^(x,t r ), i/ J (a;,i I .)G' B _i(^(a;,t r ),a;,i,.)) S /C B , j = 1, . . . ,g. 
The time-dependent Baker- Akhiezer function 

*(^».m,.)=(|;£: : :;! :: ;: :: >) ^ 

is defined by 

if)x{P,x,x Q ,t g ,to, g ) = cxp ( - -(u(a;,t fl ) - u(x ,t g )) (3.85) 
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+ i dx' <f>(P,x',t g ) - ds(-F g (z,x ,s)(/)(P,x ,s) + G g -i(z,xo,s))j. 

JXO "'to, 8 Z / 

ifj 2 (P,x,x Q ,t g ,t Q . B ) = -ipi(P,x,X ,t g ,to lS )(p(P,X,t B ), 

Pe^MPoo}, (x,x ,t g ,t , g ) 6M 4 . 
Combining relations ( 3.67a ), ( [3. 67b ), and ( 3.68a ) one concludes 



,(x, t g ) = tin ( (-l)^- 1 Mj (sf, i B ) J = -i In I (-1) B /T 1 "ifo *g) 



3=1 



3 = 1 



We will also need the following asymptotic high-energy expansion 



3=0 



4. Symmetric functions 

Let g £ N be fixed and define 

S k = {I = (4, • • • , 4) e N fe | < • • • < 4 < &}, k < 
4 J) = {&= (4, ■• • .4) e S fc | 4* ^ j}, < a - l. 



Define 



ges k 



4> 



= 1, <^V) = (-1)" E Wi-W*. * < fl - i, ^V) = o, 

where /i = (/xi, . . . , /x B ) £ C B . One recognizes the simple pattern, 

*i(m) = -2Jw, * 2< ^ = £ etc -' 



£=1 

£1 



^1,42 = 1 

<1«2 



=-£m£, $ 2 J V) = £ M4^ 2 , etc. 



=1 



£l,£2=l 



Let i? , ■ ■ • , En be iV + 1 complex numbers, where N — 2jj or N — 2g + 1 depending 
underlying hierarchy of soliton equations. For brevity we introduce 

E={E ,...,E N ). 

We will need the following elementary result. 

Lemma 4.1. For zeC such that \z\> max{|i?o|, ■ ■ ■ , |£W|} we have 



N 



n a--)) =e^ 



-1/2 



fc=0 



Co© = 1, Cfc (£) = 



3o,—,jjV=0 
30 H h3iV=fe 



Proof. It suffices to apply the binomial expansion. 

6 (2n - 1)!! = 1 ■ 3 • • ■ (2n - 1), and (-1)!! = 1. 
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The few first terms read 

N N N 

Ci(E) = - £ Ej, c 2 (E) = -J2 E 3 E k + ^ E l etc ' ( 4J ) 

3=0 j,k=0 3=0 

3^k 

Next, assuming \ij ^ fj,ji for j ^ j', introduce the g x g matrix U B (fi) by 

LMaO = 1, U^)=(—-f -) . (4.8) 

Lemma 4.2. Suppose fij ^ fij> for j ^= j' . Then 
Proof. Observe that we may write 

Using Lagrange's interpolation result, Theorem |A.1| (replacing fc by g — fc in ( |A.3j )), proves the 
result. □ 

Of crucial importance for our approach is the fact that we are able to express fj and F r in 
terms of elementary symmetric functions of fx±, . . . , fi B . The expression is given below for the 
homogeneous case only, denoted by fj and F r , where the integration constants q for I £ N vanish. 
We start with fj. 

Lemma 4.3. Let Cj(E_) be defined as in ( |4.6| ). Then we infer the following results for the KdV 
and the Toda hierarchie^\ 

j'Ag 

£=I>-*(£)**(m)- (4-11) 

For the AKNS hierarchy we obtain 

jAfl jAfl 

h = -i?X>_ fc (£)* fc (g), ftj = ij>X>-k(D*fc(K), (4.12) 

k=Q k=0 

where v_ — {v\, ... ,v s ). 

In the sGmKdV case we have 

3 3 

/; = X>-k(I)*fcQf)> = X>_fc(£)* fc (i£), i = o,...,fl-i,flGN, (4.13) 

fe=0 fe=0 

and 

/ fl = * B (M), ft s = * B (£), fl£N . (4.14) 

Proof. The proof is identical in all cases, and is based on the high-energy expansion of the Green's 
function of the corresponding linear operator L or M in the Lax pair or zero curvature formulation 
of the hierarchy considered. We provide the details for the KdV hierarchy only. 
Using Lemma |4.l| we find 

F e (z) _Y] S ]=1 {z-H) _ 1 lX B = i(l-f) 



j=0 m=0 m=0 k=0 



7 n Am = min{n, m}. 
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Combining this result with the high-energy expansion (3.13) proves the result. 



□ 



Remark 4.4. Observe that the right-hand side of (4.11) is defined for all x G M, but when we 
sample it at integer values x = n G Z, it coincides with fj(n) for the Toda lattice. Thus we have 
in some sense a continuous extension of the Toda hierarchy (cf. also Lemma 5.' 

Theorem 4.5. Let r G Nq. For both the KdV and the Tl case one derive$\ 



p=(r- fl )V0 

For the AKNS hierarchy one infers 

r r 

F r { i x j ) = -iq Y, c p {E)^l p { t ), H r { Vj )=ip J2 

p=(r-g)V0 p=(r-g)V0 

For the sGmKdV hierarchy one conclude^ 



p=( r _l_ £| )V0 
r — l 

p = ( r _l_ fl )V0 



Proof. It suffices to consider the KdV and sG cases. By definition 

r r (> — ^)Ag 

F r (z) = ^ f r -(Z l = E^ E *m(if) ( T-<-m©- 



Consider first the case r < g. Then 



r—p 



F r (z) = E C P {E) E 2^*r-^-p(M) 
p=0 



£=0 



and hence 



^(Mi) = E c f® $ - p (M) 



p=0 



using ( |03| ). In the case when r > q + 1 we find 



i? r (z) = E *mCE) E 



m=0 



p=0 

s 



= E^(£)(Ew* fl ~V~ fl ~ p + E cp(^)E^Qf)^ r - p - 

p=0 ^=0 p=r-g+l ^=0 

r— fl r r—p 

= F g (z)J2cv(E)z r - 3 - p + E Cp^E^^)^"^ 



p=0 
r—g 



p=t — g+1 



£=0 
r—p 



F B (z)Yc P mz r - B - p + E c p (£) E *r-p- m (M)< 

p— p— r— 0+1 m— 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



(4.21) 



8 n V m = maxjra, m}. 

9 Since r is independent of g, one ob tains f T = ae"'", h r = f3e lu with o,/3 £ C independent of a, /3, and 
f q , h q , q = 1, . . . ,r — I constructed as in (3.75 ). 
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Hence 

r 

F r ( H )= J2 c p(eH 3 1 p (h), (4.22) 

p=r— a+1 

by using ( A.13| ) again. 

In the sGmKdV case we first observe the identity 

F r {z) = zF r _i(.z) + / r , (4.23) 

which implies 

^ = F r ^) + t= £ --*S2i00. (4-24) 

it, u, * — ' F — a a — 

^ J nJ p=(r-l-fl)V0 



using / r = ae 411 and the trace relation (3.87). 



□ 



5. DUBROVIN EQUATIONS AND LINEARIZED FLOWS 



Dubrovin [EM made the fundamental observation that the Dirichlet divisors for the KdV equa- 
tion satisfy a first-order system of differential equations. Solving this system can then be used to 
recover the function V by appealing to a trace formula (cf. ( |5.28 )). 

Before we state the Dubrovin equations we need some notation. Let g £ N. We start by 
constructing the hierarchies as explained in Section |^. In particular, we construct the function 
F g with its zeros \i = . . . , (j, g ), and define the corresponding hyperelliptic curve K, g . (In the 
AKNS case we also construct the function H g .) Next, fix an r G No, and construct the function 
F r . The integration constants in the definition of F r are assumed to be independent of those used 
to construct F g , and to emphasize this fact we denote it by F r and the corresponding constants 
by cj>. The Dubrovin equations give the evolution of fi = (jjti, . . . , ^ g ) in terms of the deformation 
(time) parameter t r according to the rth equation in the hierarchy considered. 

The KdV hierarchy. In our setting the Dubrovin equations for the KdV hierarchy read [[l6) , 
@, Sect. 12.3, ||, ||, @, Chs. 10, 12, @, Ch. 4, 

d d 
—fj,j(x,t r ) = F r (nj(x,t r ),x,t r )— Hj(x,t r ) (5.1b) 

-_2z vSh^Ml p (u .( x t) xt) 

lU=i 3 W( X ^r) - Hl\X, t r )) 

for j — 1, ... ,g. The initial data for ( |5~l| ) on /C fl equal 

Kx,t , r ) = A (0) W, (5-2) 

where p, = (/ti, . . . , /t g ) denotes 

1 

AiOMr) = (fij(x,t r ),--F g , x (nj(x),x,t r )) e fC B , j = l,...,g. (5.3) 

We remark that (5.1a) is an immediate consequence of (f37|) and ( |U8| ), while ( pTb| ) follows from 
(pT8| ) and 

F Sl t r = F r F 3 ^ x — F ryX F g (5-4) 

upon taking z — fj,j(x,t r )- 

The AKNS hierarchy. In this case the Dubrovin equations for fi = (/2i, . . . , fi g ) are given by 

^(x,t r ) = -2z y(Aj(Mr)) (5 5a) 

dx IUjtjW{x>t r ) - fie(x,t r )) 
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d , s _ F r (nj(x,t r ),x,t r ) 9 
dt r mX '^ } ~ iq(x,t r ) dx ^ X ' tr > 

_ y(^j(x,U)) 



(5.5b) 



tr ) I1L,- (Hj(x,t r ) - fli(x,t r )) 



F r (/ij (x, £ r ) , x, £ r ) , 



for j = 1, . . . , g, with initial data on /C B 

fi(x,to,r) = P^\x), 

where 

p,j(x,t r ) = (nj(x,t r ),G e+ i(pj(x),x,t r )) G /C fl , j = 1, . . 
For the corresponding evolution of P = (Pi, . . . , P g ) we have 



d 



y(pj(x,t r )) 



v (x t r ) — —2i 

d d 
— Vj(x,t r ) = H r (vj(x,t r ),x,t r )—Vj(x,t r ) 
ot r ox 



(5.6) 
(5.7) 

(5.8a) 
(5.8b) 



y(v j (x,t r )) 



for j = 1, . 



where 



p(x,t r )]\}_ Lo (v J (x,t r ) - u t (x,t r )) 
with initial data on IC g 

0(x,t . r ) = P ( °\x), 



H r (uj (x , t r ) , x, t r ) , 



(5.9) 



(5.10) 



Pj(x,t r ) = (Uj(x,tr), -G a +l(Uj(x),X,tr)) & ^fl , j = 1, 

The Toda hierarchy. Here the Dubrovin equations for fx = (jli, ...,fi g ) read pH , | Q , [ f5l| , 
§, Ch.4, 



A M .(„,t r ) = 2 yJMhhA 

dt r 3 ' nf^j(Mj(".*r) - fJ>e(n,t r )) 



F r (fj,j(n,t r ),n,t r ), 



for j = 1, . . . , g, with initial data on /C B 



where 



A(Mo,r) = A (0) ( ri )> 
p,j(n,t r ) = (fij(n,t r ), -G 3+1 (pj(n),n,t r )) e K s , j = 1,. 



(5.11) 



(5.12) 



(5.13) 



We note that ( 5.5a ) and ( 5.8a ) form ally coincide with ( 5.1a ). The case of the Toda hierarchy, 
however, is quite different since (5.11) concerns the independence of pb(n,t r ) and no analogous 
first-order nonlinear difference equation concerning the n-dependence of fi(n,t r ) (i.e., an analog 
of ( 5.1a| ) or (5.5a)) appears to be known. In this context we refer the reader to Lemma 5.6, where 
we continue this discussion. 

The sGmKdV hierarchy. Finally, in the case of the sGmKdV hierarchy the equations for n 
read §1 



y(fij(x,t r )) 



Hj,t r (x,t r ) = 2 



nLlOjOMr) -ln(x,t r )Y 

F r (nj(x,t r ),x,t r ) y(fjj(x,t r )) 

Hj(x,t r ) Y\ 3 i = i{^j{x;t r ) - fj,t(x,t r )) : 

j = 1, . . . ,0, (x,t r ) e Q, 



(5.14) 
(5.15) 



with initial data 



fij(xo,to,r) EfC g , j = 1, . 



(5.16) 
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The corresponding equations for v_ equal 



( , \ _ 9 gr^(MrQ 1 Mr) V^j{x,t r )) 



j = 1, . . . ,0, (a;,i r ) e 



with initial conditions 



%(#o,*o,r) e £„, j = l,...,g. (5.19) 

Next we will prove that the Abel map provides a clever change of coordinates that linearizes the 
Dubrovin flows. This will turn out to be a consequence of the fact that F r ([ij) can be expressed 
as a linear combination of the functions . Using Theorem |4.5| it is immediate that this is not 
only the case for the KdV hierarchy, but also for all the other hierarchies discussed in this paper. 

Theorem 5.1. Suppose /j,(x,t r ) = (/Xi(ic,t r ), . . . ,(i g (x,t r )) satisfies the Dubrovin equations (5.1) 
and assume that fij ^ (Xji for j ^ j' . Let r G No and introduce 

rAg 

F r (fX j )=J2 d r,k®% ) {tJ L ), dr,0,... > dr,rA S eC. (5.20) 
fc=0 

Then the Abel map 

Ap {P>j{x,tr)) = ^^(Aj (Mr)), •••,4i^, fl (AK a; > **-))) (5.21) 

linearizes the Dubrovin flows (|5.1|) in the sense that 



id hencef®\ 



g S fl 

gf^Apbtifafatr)) = ~ 2i c k,£d r , s -e (5.22) 

r 3=1 «=iv(n-r) 



QLp (Pii{x,tr)) = OL Po {V^ Xo ^ to r) ) - 2i(x - x )c k , s d fi - 2i(t r - t . r ) c k ,td riB -£ . (5.23) 

fc=lV(fl-r) 

Proof. One computes, 

^ S g a rH(x,t r ) 



ot r j=1 cn r . =i j Po 



^ ^ JPn V{P) 

r 3=1 t=X J ^> »\ / 

__9-V~ v V~ v f-j(x,t r ) e 1 y(fij(x 1 t r )) ~ 

" j^[t=X (M^) rim/j (Mj <r) - Mm (z, tr )) 

fl fl B 

= -2i} j }^Ck,eU B (n(x,t r ))ijF r (nj(x,t r )) = -2i } c k ,id r:S ~£, (5.24) 

3=1 <=i <=lV(a-r) 



using Lemma 4.2 in the final step. As for the ^-variation, we observe that in-derivative of fXj 
coincides with the x-derivative in (5.1), and hence it is a special case of (5.24). This proves the 
theorem. □ 



10 The situation here resembles the one in classical mechanics where, by a canonical change to cyclic coordinates, 
the momentum pj becomes a constant of motion and thus qj(t) = <Jj(to) +Pj(t — *o) is linear in time. 
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Corollary 5.2. The Abel map linearizes the Dubrovin flows for the KdV, AKNS, Tl, as well as 
the sGmKdV hierarchies. 



Proof. Th eorem 4.5 shows t hat F r {jij) (and H r (vj) in the AK NS c ase) indeed satisfies the as- 



sumption ( 5.20 ) of Theorem 5T, and hence the key calculation ( 5.24 ) carries over to the AKNS, 
Tl, and sGmKdV systems. The special case r = gives the ^-variation in all but the sGmKdV 
case which, however can easily be verified by explicit computation. □ 



Remark 5.3. We provide a few more details in the AKNS case. Suppose fi{x, t r ) satisfies ( |5.5| ) and 
similarly, v_(x,t r ) satifies (5.8), with jij ^ jiji and Vj ^ vy for j =/= j'. Let r € No and introduce 

rAg rAg 

F r { H ) = -iq^dr^Hv), H r {vj) = ipJ2 e r,k^\E)- (5.25) 

k=0 k=0 

Then ( [5^22] ) and ( |5~23| ) hold. In addition, one obtains the following results for the analog of 
Neumann divisors v(x, t r ). 

d 8 s 

^2J--Po,fc(^'( x >*'-)) = ~ 2i c k,£e r , B -e (5.26) 

r 3=1 fclV(g-r) 

and hence 

8 

&p {' D o(x,t T )) = ap (%i i„ !r )) ~ 2i(x - x )c kiB e . - 2i(t r - t . r ) c k je r , s -e- (5.27) 

fclV(g-r) 

Necessary and sufficient conditions on Lax pairs to linearize the flow t — > L t on J(C), where {L t } 
represents a dymamical system on the Jacobi variety J(C), with C the underlying spectral curve, 
have been considered by Griffiths j3^]]. While he considers Lax equations within a cohomological 
framework, our approach is much more modest in scope but in turn reduces the linearization 
problem to an elementary exercise in symmetric functions. 

Solving these equations we can recover the solution of the integrable equation using trace 
formulas. For the KdV hierarchy we have the classical trace formula 

20 8 

V(x,t r )= ^£ ro -2^Mj(Mr). (5.28) 

m— j—l 



For the AKNS hierarchy we have 

Px (^5 ) 
p(x,tr) 

q x (x, t r ) 
q(x,t r ) 



2g+l b 

i E m - 2iy^ j fj(x,t r ), (5.29a) 

rn—O j — l 

28 + 1 g 

-i ^2 E m +2i^jjL j (x,t r ), (5.29b) 



m=0 j=l 

while for the Toda hierarchy one obtains 



t 2 = _1 G B+ i(/jj(rt,t r ),n,t r ) 

2 ~i Uk&ih ■("•>*r) - Hk{n,t r )) 



-l b (n,t r ) 2 -^J2f, 3 (n,t r ) 2 + ^J2 E ^ ( 5 - 30a ) 

j — 1 m— 

1 2g+l g 

6(n,*r) = -- 2 ^m+X^M'-)- ( 5 - 30b ) 



2 

m— j—l 
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The "trace" relation for the sGmKdV hierarchy (it would be more appropriate to call this a 
"determinant" relation) was given in ( |3.87| ), 

u(x, t fl )=i In ((-lJ'a^n^M,)) = -iln ^(-l)^" 1 JJ ^(^^)) ■ (5-31) 



3=1 



Remark 5.4. It is important to observe that if one postulates the Dubrovin equations (5.1), and 



defines V using the trace formula (5.28|), one could show by a long and tedious calculation that 



V indeed satisfies the rth KdV equation with the correct initial condition. The same applies, of 
course, to the AKNS, Toda, and sGmKdV hierarchies. 



Remark 5.5. For simplicity we assumed (j,j(x,t r ) ^ fiji(x,t r ) for j =/= j' in Theorem 5.1. In 
the self-adjoint cases, where {-E m }m=o,...,JV C R, this condition is automatically fulfilled for all 
(x,t r ) £ R 2 since then all fij(x,t r ) are separated from each other by spectral gaps of L(t r ) or 
M(t r ). In the general nonself-adjoint case this is no longer true and collisions between the fij's 
become possible. Nevertheless the Dubrovin equations, properly desingularized near such collision 
points, stay meaningful as dem onstr ated in deta il by Birnir 0, jl0| in the case of complex- valued 
KdV solutions. In particular, (5.23) (and (5.27)) remain valid in the presence of such collisions 
due to the continuity of a Po (-). 

We already mentioned in the paragraph following ( pM3| ) that the Dubrovin equations for the 
Toda hierarchy differ from the ones associated with the remaining soliton hierarchies in the sense 
that they do not seem to govern the n-dependence of fi(n, t r ). We now show how to use Theorem 
5.1 to obtain a first-order Dubrovin system for p,(x, t r ) in x, whose solution coincides with fi(n, t r ) 
at the integer points x = n S Z. Since the independence of (1 plays no role for this argument, we 
ignore this dependence in the following result. 



Lemma 5.6. Abbreviate A — (A\ 1 . . . ,A g ) = A P ^ 
( 3.()2lf ). Consider the Dubrovin-type system 



(Poo + ) and £t(no) — p,^, with fi(n) defined in 



= 1 



/t(n ) = A j 
with c m . n defined in ( pjj| 
ti Q {x) at integer values x 



(5.32a) 
(5.32b) 

Denote the solution of ( |5.32| ) by A ( a; )- Then fj(n) coincides with 
n G Z, that is, 



/t(n) = A ( n )j n G 2- 
Proof. First we recall the well-known result (see, e.g., pj, Sect. 3) 

~ "p (^A(no)) = ( n ~ n o)A Poo _ (^00+) = (n - n )A. 
In order to complete the proof we only need to stablish that fi Q (x) satisfies 

QLpoi^^x)) -a Po (V^ Xo) ) = (x - x Q )A Pcxi _(P oc+ ) = (x-xq)A 
that is, we need to show 

A. 



(5.33) 



(5.34) 



(5.35) 



(5.36) 



But equation ( 5.36| ) follow s immediately from ( 5.11) ( identifying t r with x and ignoring its n- 
dependence), (^.20[) , (|5~2^ ) (multiplied by i), and (|5.32|). □ 



Thus, the solution p- Q (x) of (5.32) provides a continuous interpolation for jl(n). In fact, it was 
our attempt to prove a result like Lemma 5.6 which led us to reconsider Dubrovin equations and 
ultimately resulted in Theorem 



4.5 



the explicit connection between $jr^ (p) and the polynomials 
) defining the hierarchy in question. 
Toda systems as integrable discretizations of continuous systems are also studied in H . 
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6. THETA FUNCTION REPRESENTATIONS 

The fundamental problem in the construction of algebro-geometric solutions of soliton hierar- 
chies is the following. Fix g G N, and pick a stationary solution of the gth equation in the hierarchy. 
Then solve the rth time-dependent equation (for any r S N) with the given stationary solution 
as initial datum and express the solution in terms of the Riemann theta function associated with 
IC B . In the KdV case this procedure has been pioneered by Its and Matveev in their celebrated 
paper pp| . Subsequently, the algebro-geometric approach to integrable equations was developed 
in papes by Date, Dubrovin, Krichever, Matveev, Novikov, Tanaka, and others fig ], ]ri| , p^ j, [ p3| , 
p4j , p6| ], Q, p3| , Q. We briefly recall the results for the hierarchies studied in this paper. 
Detailed discussions for the KdV, AKNS, and Toda hierarchies as well as the sG equation and 
other completely integrable systems can be found, for instance, i n m , Chs. 3, 4, 6, Oh Sect. 
12.4, H, H, ||, @, @, H, f§, @, @, Sect. 4.4, @, jp, [||, @, @|, ill, lH, Ch. 
4 and the references therein. (Without explicitly repeating this in each case below, we exclude 
collision points for fj,j(x,t r ), that is, we will always assume fij fj,ji for j ^= j'.) 

The KdV hierarchy. Let V^(x) be a stationary solution of the gth KdV equation associated 
with divisor T>^o)r x y We seek the rth KdV flow 

KdV r (V)=0, V(x,t , r ) = V {0) {x), xeR. (6.1) 
The solution is given by the Its-Matveev formula ]4C| | 

V(x,t r ) = A-2dl\n(0(z(fi(x 1 t r )), (i,( r )el 2 , (6.2) 

where 

A(as,*o,r) =A (0) (^) (6-3) 

and 

z([i{x,t r )) = a Pb (X»^( X)tr) ) - A Po (Poo) + S Po , (6.4) 

with Sp the vector of Riemann constants and A a /C B -dependent constant. 

The AKNS hierarchy. Let (p^(x), q(°\x)) be a stationary solution of the gth AKNS equa- 
tion associated with divisors T> ^o)/ x \ and T>.{ )^ x y We want to construct the rth AKNS flow 

AKNS r (p,q) = 0, (p(x,t ,r),q(x,t , r )) = (jpW(x),qW(x)), x e M. (6.5) 
The solution reads 

/ , x / , ^(z_(^(x ,^ r ))) fl(z + (j>(s, £,.))) 

p(a;, t r ) = p{x , t ,r) T , — Try — ttt^ — 7T7 - ^ x (6.6a) 

x exp(— 2i(a; — x )eo — 2i(£ r — t a r )e r ), 

Q(X,t r ) = g(x ,i ,r)^7 7^7 — ttt^t — 7^7 — rrrr x (6.6b) 

x exp(2i(x — a; )eo + 2i(i r — to r )e r ), 

(z,i r ) 6 R 2 , 

where e r is a /C fl -dependent constant, and 

^(i,to, r ) = r(i), £0Mo,r) = ^ (0) (aO (6.7) 

and 

2 ± (Q) = ap (X'Q)-4po( i3 oo ± )+Hp , Q=(Q 1 ,...,Q,), (6.8) 

with eo, e r are ^-dependent constants. 

The Toda hierarchy. Let (a^°\ b^) be a stationary solution of the gth Toda equation 
associated with divisor £>^<o) / n \ . We are interested in the solution of the Tl r flow 

Tl r (o,6) = 0, (a(n,t Q>I .) > 6(»,tQ, I .)) = (a (0) (n), 6 (0) (n)), n € Z. (6.9) 
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Here the solution equals 

(z,(A(n - 1, t r )))0{z + {jl{n + 1, tr))) \ 1/2 



a(n, t r ) = a 



0(z + (gn,t r W 

2fl + l B 



b(n, t r ) = Aj - -= E ™ ~ c ^ s) 7w ln 



J'=l 



m=0 



i=i 



9(w + z+{Kn,t r ))) 
{w + z + (ji{n - Mr))) 

(n,t r ) G Z x R, 



w=0 



(6.10a) 
(6.10b) 



where 



and 



fi(n,t ,r) = /t (0) (n) 



z + {gn,t r )) = tt Po ^(„, ir) ) -Ap (Poo+) +S Po , 



(6.11) 



(6.12) 



with ii/0a /C g -dependent constant. 

The sGmKdV hierarchy. Let (x) be the solution of the gth stationary sGmKdV equation, 
that is, 



subject to the constraints 



9 s -i, x - i(Pe mW -ae- mlO> ) = 0, g G N, 



2g 



(6.13) 



(6.14) 



Let r G No. We are seeking the solution u of sGmKdV r (u) = with u(x, to,r) = w^(x). Here the 
solution reads 



u(x, t r ) — u(xq, t r ) 



(6.15) 



2iln 



9(z(fi(x, t r )) + A)d(z(fi(x , t , r ))) 



J{z(Mx ,t 0i r)) + A)6(z(jl(x 
where P = (0,0), A is a half-period, 

A = A Po (P oc ), 

e a /C fl -dependent constant, and 

= «Po (%) - ^Po (^oo) + H Po 

The linear equivalence of 'Dp oa p,{x,t r ) an< ^ T~ > Qov{x,t r )i that is, 

Sfb(%(x,i r )) =Sfl,(%,t r )) + A, 

shows that 

z(v(x,t r )) — z(p,(x,t r j) + A. 



-^yyexp(-ie (a;-a;o)) I , 



(6.16) 



(6.17) 



(6.18) 



(6.19) 



In each of the theta function representations of this section one should keep in mind that tha 
change of coordinates effected by the Abel map straightens out all Dirichlet and Neumann flowsE] 
on the Jacobi variety J(/C B ) of /C B , that is, z((l(x, t r )) and z{v_{x,t r )) in ( |6.2[ ), (6.6), ( |6.1C| ) and 
( |6.15 ) are linear in x (respectively, n) and t r . 
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7. Examples 

KdV and sG. Pick E = and E\,..., E 2s G C, E m ^ E m ' for m ^ m! and solve 

± N ( x ,t a ) = -2i — B - yp^Mt — (7 .i a) 

9 9 1 fl 

— /i 3 -(a!,t a ) = —fXj(x,t B ) 16 a JJ^, j = l,...,fl, (7.1b) 
with Q = rim=i and i? 2B +i(z) = z ]lm=i( z ~ #ro). Define 

£1 

uOMg) =iln(Q- 1/2 n^(^^ fl )), (7-2) 

2fl B 

V(z.* B ) = ^^ m -2^^(^* B )- (7-3) 

m— J — 1 

Then u and 1/ satisfy the sG equation and gth KdV equation, respectively, that is, 

41^,4, = sin(u), KdV s (F)=0 (7.4) 

for the following choice of q, 

5o = 1 ' 5l = ( 16Qi/2 ~ Cl( ^' c e = -J2<h>ce-p(E), e = 2,...,g. (7.5) 

The isomorphism between algebro-geometric KdV and sG equations is of course well-known and 
has been discussed, for instance, in ||, Q. 

AKNS and Tl. Pick E Q , . . . , E 2s +i G C, E m ^ E m , for to ^ to' and solve 

£*A*,tr) = -« n . ff^» (7.6a) 

—flj(x,t r ) = —^fij(x,t r ) d n®nHp), j = l,...,S (7.6b) 

n=(r-fl)V0 

— Vj(x,t r ) = — 2i = — - — — ^ — — , (7.7a) 

c>x n n^i(^(^*r)-^(^,M) 

■g-Vj(x,t r ) = -— Vj {x,tr) Yl i = (7.7b) 

n=(r-fl)V0 



and 



Px(x, t r ) 



where i?2 B +2(z) = Ilm=o ( 2 _ ^m) and rf „ G C. Define 

20 + 1 fl 



p(x,t r ) 
q(x,t r ) 



i Y E m -2iJ2"j( x ^r), (7.8a) 

m— J — 1 

2fl+l £1 

-i ff m + 2i^/j, j (x,t r ) (7.8b) 

m— i=l 



and 



a(n, ir ) 2 =^ v(Ai(n.*r)) 



E[fc^j(Mj(«>*r) -fJ.k{n,t r )) 
1 1 fl i 2b+1 

6(n,i r ) 2 --^ Mj (n,t r ) 2 + - ^ £*, (7.9a) 



4 

^' — 1 m— 
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1 2g+l g 

b(n,t r ) = - j £ ^ m + ^^(n.*r). (7.9b) 

m=0 j — 1 

Then (p, g) and (a, 6) satisfy the rth AKNS equation and the rth Toda lattice (Tl) equation, 
respectively, that is, 

AKNS r (p, q) = 0, Tl r (a, b) = (7.10) 



for the same choice of eg, I = 1, . . . , r in both equations (7.10) (depending on the choice of d n in 
(pb|), (|7Jbl)). 



Remark 7.1. These examples provide interesting connections between the KdV B and sG equation 
(where N is even and S {£ , m } m =o,...,A r )j an d AKNS and Toda hierarchies (where A" is odd), 
respectively, and illustrate the fundamental role of the Dubrovin equations as the common un- 
derlying principle for hierarchies of soliton equations. In particular, our approach establishes an 
isomorphism between the classes of algebro-geometric solutions of these pairs of integrable sys- 
tems. Indeed, once the hyperelliptic curve K, g is fixed, algebro-geometric solutions of the KdV B 
and sG (respectively, algebro-geometric solutions of the rth AKNS and rth Tl equation) are 
just certain symmetric functions (i.e., "trace" relations) of the solutions ni{x, t g ),..., ^ g {x, t g ) 
(resp. ni(x, t r ), . . . , /i g (x, t r ), v\(x, t r ), . . . , v B (x, t r )) of the corresponding Dubrovin equations on 
JC g . Analagous considerations apply to the nonlinear Schrodinger equation and the (continuum) 
Heisenberg chain (see, e.g., |f27j| ). The interesting problem which types of symmetric functions of 
/Lii, . . . ,/i B (i.e., which types of "trace" formulas) actually lead to completely integrable hierarchies 
is currently under investigation. 

Appendix A. Lagrange interpolation formulas 

In the following we suppress the (x, ^-dependence as it will be of no importance in this ap- 
pendix. 

Fix q £ N and recall that 

s 

F 3 (z) = l[(z-» j ), (A.l) 

i=i 

which implies that (F^ = dF g /dz) 

e 

F'M) = ]J^k-^). (A.2) 

3=1 

The general form of Lagrange's interpolation theorem then reads as follows. (For convenience of 
the reader we supply its proof even though the result is well-known.) 

Theorem A.l. Assume that fit, . . . , /i g are g distinct complex numbers. Then 

m=l,...,fl+l, fe = 0,...,s-l. (A.3) 

Proof. Let Cr be a circle with center at the origin and radius R that contains the zeros /ij of the 
polynomial F g and which is oriented clockwise. Cauchy's theorem then yields 

i r /-m—1 ^m— 1 ,,m—X 

— i dc - = - + V ^ 

2mf CR ^F g (0(C-z) F g (z) Z^F'^iM-zy 

z± /xi,...,/i g! m= 1, ...,0 + 1. (A.4) 

However, by letting R — > oo we infer that 

1 f Am — 1 R m ~^ 

T~ f rf C rp , t \(/- \"i im FTm m = l,...,fl + l, (A.5) 

2 ™ Jc R F g (C) (C Z) R-ooF, (R) 
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which implies 



171 



Using the symmetric functions tyj we may write 



= F B (z)6 mt3+1 . 



and 



3=0 



z-Hj 



(A.6) 



(A.7) 



(A- 



fe=0 



Expanding both sides of equation (A.6) in powers in z, using (A.7) on the right-hand side and 
flA.8|) on the left-hand side, proves (A.3). □ 



The simplest Lagrange interpolation formula reads in the case k = 0, 

• 0- 



B m-1 

> T = 5 m „, 171=1, 



(A.9) 



Fo r use in the main text we finally observe the following equalities. Adding (A.7) to fj,j times 
(O) we find 



F 3 (z) + N = £ 1 + rf ) 



B-l 



k=0 



However, we also have 



F (z)+(Xj 



s-i 



fe=0 



using QA.SD and recalling = 0. Thus we conclude 

+ Mi^Qf) = *Si(£f), * = 0, • ■ ■ ,3 - 1. 

Finally, we will show 

fc 

^ Al ^ fe _,(/£) = $i i) (M), fc-O,...,0, 
by induction. Equation ( A.13Q clearly holds for k = 0; next assume that 

holds. Then 

fe k 



(A.10) 



(AH) 



(A.12) 



(A.13) 



(A.14) 



(A.15) 



e=o 



fc-i 



using first the induction hypothesis and then ( A.12j ). 
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